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ABSTRACT 

A model for heliospheric solar wind charge exchange (SWCX) X-ray emission is applied to 
a series of XMM-Newton observations of the interplanetary focusing cone of interstellar helium. 
The X-ray data are from three coupled observations of the South Ecliptic Pole (SEP, to observe 
the cone) and the Hubble Deep Field-North (HDFN, to monitor global variations of the SWCX 
emission due to variations in the solar wind) from the period 24 November to 15 December 
2003. There is good qualitative agreement between the model predictions and the data with the 
maximum SWCX flux observed at an ecliptic longitude of ~ 72°, consistent with the central 
longitude of the He cone. 

We observe a total excess of 2.1 ± 1.3 LU in the O VII line and 2.0 ±0.9 LU in the O VIII line. 
However, the SWCX emission model, which was adjusted for solar wind conditions appropriate 
for late 2003, predicts an excess from the He cone of only 0.5 LU and 0.2 LU, respectively, in 
the O VII and O VIII lines. We discuss the model to data comparison and provide possible 
explanations for the discrepancies. 

We also qualitatively reexamine our SWCX model predictions in the ^ keV band with data 
from the ROSAT All-Sky Survey towards the North and South Ecliptic Poles, when the He cone 
was probably first detected in soft X-rays. 

Subject headings: X-rays: observations — Solar System 


1. Introduction 

Solar wind charge exchange (SWCX) X-ray 
emission in the heliosphere originates when a high 
charge state ion of the solar wind interacts with a 
neutral atom and gains an electron in a highly ex- 
cited state which then decays by emission of an X- 
ray with a characteristic energy of the ion. SWCX 
emission occurs over a wide range of spatial scales 
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including comets (where such emission was first 
conclusively observed, e.g., Lisse et al. 1996), the 
near-Earth environment where the neutrals are ex- 
ospheric material in and near the magnetosheath 
(e.g., Cravens et al. 2001; Snowden et al. 2009), 
and throughout the solar system where the target 
atoms are from the interstellar medium (e.g., Cox 
1998; Cravens 2000; Smith et al. 2005). SWCX 
emission comprises a significant part of the diffuse 
X-ray background observed at energies < 1 keV, 
and can be both a signal of interest in studies of 
the solar system or an unfortunate contamination 
component in studies of the hot, diffuse plasmas 
in the Milky Way and beyond. 

In theory. SWCX emission exhibits very dis- 
tinct spectral characteristics with respect to ther- 
mal plasma emission. For instance, the line inten- 
sity ratios in triplets of the Hc-like ions O VII and 
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Ne IX are very different and could be used as spec- 
tral diagnostics to separate the two mechanisms 
(e.g., Lallement 2008, and references therein). Fu- 
ture missions which include non-dispersive mi- 
crocalorimctcr spectrometers will be able to sep- 
arate thermal from SWCX emission using the 
triplets (e.g., Snowden 2008). However for dif- 
fuse sources, with the spectral resolution of current 
missions it is impossible, in general, to distinguish 
individual SWCX X-rays from those with a more 
distant origin, and alternative diagnostic informa- 
tion must be used to separate the components of 
the X-ray background. Temporal variation can be 
used to a certain extent depending on the length of 
an observation, and indeed this method was used 
in cleaning the ROSAT All-Sky Survey (RASS, 
Snowden et al. 1997). SWCX emission from the 
near- Earth environment will track relatively well 
variations in the solar wind flux as measured by 
the Advanced Composition Explorer (ACE. Stone 
et al. 1998) and WIND (Ogilvie et al. 1995) in- 
struments. However, identifying and modeling a 
temporally variable component can only address 
a fraction of the SWCX emission. Emission from 
interactions with interstellar (IS) neutrals drifting 
through the solar system will typically have mini- 
mal temporal variation as the observed flux along 
any line of sight is an integral over weeks of so- 
lar wind conditions. Because of the long (relative 
to a typical X-ray observation) times of any ex- 
pected variation, IS SWCX effectively contributes 
a zero-level offset to observations of the diffuse X- 
ray background of more distant origin. 

While there are models for SWCX emission 
from the heliosphere, testing them is problematic. 
Koutroumpa et al. (2007), using a self-consistent 
model of the heliospheric SWCX emission, man- 
aged to associate observed discrepancies in XMM- 
Newton and Suzaku observations separated by sev- 
eral years with solar cycle-scale variations. A re- 
cent paper (Snowden et al. 2009) was relatively 
successful at using an XMM-Newton observation 
to test a model (Robertson & Cravens 2003) for 
emission from the near- Earth environment. In this 
paper we use a different series of XMM-Newton 
observations to search for emission from the solar 
system. Specifically, we use multiple observations 
of the south ecliptic pole as Earth passes over the 
the ISM helium focusing cone downstream of the 
Sun to search for a correlation with the expected 


variation from the model. 

2. Observations 

2.1. Helium Focusing Cone 

The flow of interstellar neutrals through the so- 
lar system is due to the motion of the heliosphere 
at about 25 km s~ 1 through the local IS cloud 
(LIC). This material, a gas of mostly hydrogen 
atoms with about 15% helium flows from the direc- 
tion of A, (3 ~ 252°, 9° (o, 6 ~ 252°, -14°). This 
places the Earth upstream of the Sun in the inter- 
stellar neutral flow in early June and downstream 
in early December every year (Gruntman 1994). 
Although both radiation pressure and gravity af- 
fect the hydrogen trajectories (e.g., Quemerais 
et al. 1999; Lallement 1999), only gravity sig- 
nificantly affects the helium trajectories which 
execute Keplarian orbits and form a “focusing 
cone” downstream of the sun, resulting in a lo- 
calized downstream enhancement of helium ob- 
served annually by Earth orbiting and LI space- 
craft (Bzowski et al. 1996; Frisch 2000). 

Many diverse sets of observations have con- 
firmed this helium focusing cone phenomenon, and 
have put strict constraints on the helium flow pa- 
rameters (density, temperature, velocity vector), 
including direct neutral gas measurements (Witte 
2004) and UV backscattering (Lallement et al. 
2004; Vallerga et al. 2004) as well as pickup ions 
(Gloeckler et al. 2004) which form when the neu- 
tral helium becomes ionized primarily from solar 
photons and subsequently “picked-up” by the so- 
lar wind flow (Moebius et al. 1985). Of these, 
pickup ions have been particularly useful because 
they display a clear spectral signature, a cut-off 
at two times the solar wind speed resulting from 
their low initial speed (~25 km/s) in comparison 
to the solar wind speed (~400 km/s). However, 
their low count rates also necessitate aggressive 
averaging, commonly 30-day running averages, to 
smooth out fluctuations. 

Recently. Collier et al. (2004), based on an 
analysis of several diverse data sets, proposed the 
presence of a secondary stream of neutral atoms 
at 1 AU located at a higher ecliptic longitude 
by 10 to 40° than the primary stream described 
above (see also Wurz et al. 2004). These neutrals 
seem to be far more energetic than the nominal 
He flow. Very recent observations from the pair of 
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STEREO spacecraft (Wang et al. 2008) showing 
two peaks in the 4-20 keV neutral hydrogen atom 
flux, one primary peak at 245° close to the nominal 
upstream direction and another secondary peak 
shifted 16° toward higher ecliptic longitudes than 
the nominal upstream direction, appear to con- 
firm this prediction. Additionally, Mars Express 
has observed a neutral atom signal consistent with 
this secondary stream proposal (Holmstrom et al. 
2008). Nevertheless, all observations of slow neu- 
trals show only one cone, and the secondary cone, 
should it exist, must consist of high velocity neu- 
trals which are not well focused. 

Because high charge state solar wind ions will 
charge exchange with the enhanced helium in the 
downstream region emitting soft X-rays, the he- 
lium focusing cone should be observable in soft 
X-rays (Cravens et al. 2001). Indeed, the ROSAT 
All-Sky Survey (RASS) data (Snowden et al. 1995) 
when observing within two degrees of the south 
ecliptic pole show enhanced emission in early De- 
cember (see §2.4) with some evidence for a sec- 
ondary peak about ten days later (or ten degrees 
higher in ecliptic longitude) than the expected lo- 
cation of the helium focusing cone (see Figure 5). 
These ROSAT data in part motivated the XMM- 
Newton observations of the helium focusing cone 
and particularly the day of year of the third ob- 
servation. 

2.2. X-ray Data 

Three pairs of XMM-Newton calibration obser- 
vations, matched pointings of the south ecliptic 
pole (SEP, to observe the focusing cone) and of 
the Hubble Deep Field - North (HDFN), were ap- 
proved by the project and scheduled for late 2003 
(see Table 1 for the observation details). The cou- 
pled observations of the HDFN were included to 
serve as controls to monitor the SWCX emission 
variation not associated with the focusing cone 
(e.g., due to variations in the solar wind flux and 
composition) which could arise from Earth’s mag- 
netosheath or more generally from the heliosphere. 

The orbit of XMM-Newton is highly elongated 
(perigee ~ 10 4 km. apogee ~ 10 5 km) and in- 
clined (inclination ~ —40°) relative to the equa- 
torial plane. In late November and early Decem- 
ber. the time period when Earth is closest to the 
focusing cone and therefore the time for these ob- 
servations. the apogee lies in the anti-solar direc- 


tion. This is fortuitous for two reasons: since ob- 
servations take place away from perigee any ob- 
served SWCX emission from the magnetosheath 
(well down the flanks of the magnetosheath and 
away from the sub-solar point) is significantly re- 
duced (e.g., Robertson & Cravens 2003) and the 
likelihood of soft proton contamination is also re- 
duced (Kuntz & Snowden 2008). 

We reduced the EPIC data using the XMM- 
Newton ESAS 1 analysis package (Snowden & 
Kuntz 2006) as demonstrated in Snowden et al. 
(2008). The data were first screened for varia- 
tions in the light curve which removed excesses in 
the normal internal particle background and more 
commonly contamination by soft protons (Kuntz 
h Snowden 2008). Next, data from the full field 
of view were extracted after the exclusion of point 
sources to a limit to 10" 14 ergs cm" 2 s" 1 . Finally 
model particle background spectra were produced 
for subtraction during the spectral fitting process. 

After extraction the spectra were then fit using 
a model which included two thermal components 
(an unabsorbed ~ 0.1 keV model for local emission 
and an absorbed ~ 0.3 keV model for the more 
distant disk and halo emission) and an absorbed 
power law representing the cosmic background, 
monochromatic lines at 1.49 keV and 1.75 keV 
representing the internal fluorescent Al Ka and 
Si Ka background, and a power law not folded 
through the instrumental response representing 
any residual soft proton contamination left after 
the screening process. Xspec APEC thermal spec- 
tral models with variable abundances were used 
where the abundance of oxygen was set to zero 
and the other abundances were fixed at 1 (based 
on the assumption that oxygen will dominate any 
SWCX emission in the XMM-Newton bandpass). 
Monochromatic lines at 0.57 keV and 0.65 keV 
were added to the model to represent the observed 
oxygen emission from both SWCX and more dis- 
tant cosmic emission, the latter which produces 
a constant contribution. The absorption of the 
halo thermal components and the power law were 
fixed at the Galactic values. Two other mono- 
chromatic lines were also added. C VI at 0.46 keV 
and another O VIII line at 0.81 keV. The addi- 
tion of the C VI line marginally improved the fit 
but left the \ 2 value unchanged with detections 


1 http://hetusarc.gsfc. nasa.gov/docs/xmm/xmmhp_xmmesas.html 
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Table 1 

XMM-Newton Observation Details and Fitted Line Fluxes in LU 
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at the < 2 a level. The O VIII line at 0.81 keV 
was insignificant. A spectrum derived from the 
ROSAT All-Sky Survey (RASS) data was simulta- 
neously fit with the oxygen emission coming from 
the thermal models rather than the monochro- 
matic lines. (The RASS data primarily constrain 
the low-temperature thermal component.) 

All of the SEP data and all of the HDFN data 
were fit simultaneously allowing only the SWCX 
oxygen Gaussians, instrumental A1 and Si Gaus- 
sians, and soft-proton power law components to 
vary between observations. The model compo- 
nents representing the cosmic background were as- 
sumed to be constant. Table 1 also lists the fit- 
ted values for the O VII and O VIII flux in line 
units (LU, photons cm -2 s _1 sr" 1 ). Also included 
are results from the spectral fits from two of the 
HDFN spectra from Snowden, Collier, k Kuntz 
(2004). The data are also plotted in Figure 1. 

Figure 2 shows the spectral fits of the SEP data 
(the fits to the HDFN data were similar). The 
fits are reasonably good for the cosmic background 
with a xt value of 1.24 for 902 degrees of freedom 
(xt value of 1.25 for 771 degrees of freedom for 
the HDFN data). Of note in the plot are the sig- 
nificantly better statistics of the PN data with the 
O VII line clearly seen at E ~ 0.57 keV and with 
the O VIII line forming the E ~ 0.65 keV hump 
on the shoulder of the O VII line. 

The fitted fluxes as plotted in Figure 1 show 
significant scatter, and any expected trend of the 
first SEP (off-cone) pointing being lower while the 
last two pointings (on-cone, the first at the nom- 
inal cone position, the second at the secondary 
cone position of Collier et al. 2004) being higher 
is completely obscured. (In fact the general trend 
is in the opposite direction.) However, the HDFN 
results which are nearly unaffected (see below) by 
the helium focusing cone enhancement, and should 
be roughly constant except for any variation in 
the solar wind flux, also show the same trend. We 
therefore used the HDFN data to normalize the 
SEP data. As noted above, the observations took 
place while the XMM-Newton satellite was on the 
opposite side of Earth from the Sun minimizing 
the SWCX emission from in and near the magne- 
tosheath (Robertson k Cravens 2003 ) allowing the 
observed SWCX X-ray emission to be dominated 
by the heliosphere. 

With this assumption we used the data from the 
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Observation 

Fig. 1. — Fitted values for the O VII (upper group- 
ing of points) and O VIII line fluxes (lower group- 
ing of points). Data points with an X symbol 
are in the direction of the SEP and data points 
with a circle symbol are in the direction of the 
HDFN. The pairs of points at observations 1—3 are 
ObsIDs (0162160101, 0162160201), (0162160301, 
0162160401), and (0162160501, 0162160601), re- 
spectively. The observation 4 data points are from 
two HDFN spectra (0111550201, 0111550401) 
from Snowden, Collier, k Kuntz (2004). 



Energy (keV) 

Fig. 2. — Spectral fit for the SEP spectra. The 
upper set of data and curves are the three PN 
spectra, the middle set of data and curves are the 
six MOS1 and MOS 2 spectra, and the bottom is 
the RASS spectrum. 
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Table 2 

Normalized 1 " 1 Line Fluxes in LU 


Obsld 

Normalized'- 1 
O VII Flux 

Normalized b 
O VIII Flux 

Normalized 0 
O VII Flux 

Normalized 0 
O VIII Flux 

0162160101 

6.5T0.8 

1.6T0.3 

7.2T0.8 

1.7T0.3 

0162160301 

8.5T1.0 

3.4T0.8 

9.3T1.0 

3.7T0.8 

0162160501 

8.6T0.7 

3.5T0.4 

8.6T0.7 

3.5T0.4 


a Data are scaled for the observed variation of the HDFN flux and then 
normalized by the model HDFN SWCX flux. 

b Scaled only by the observed variation of the HDFN flux. 

c Sealed by the observed variation of the HDFN flux and then sealed 
again by the model variation of the HDFN SWCX flux. 


Table 3 

Offset Subtracted and Normalized 3, Line Fluxes in LU 


Obsld 

Normalized 15 
O VII Flux 

Normalized 15 
O VIII Flux 

Normalized 0 
O VII Flux 

Normalized 0 

O VIII Flux 

0162160101 

2.7T0.6 

0.5+0.2 

3.0T0.6 

0.6T0.2 

0162160301 

4.5T1.0 

1.6T0.7 

4.7T1.0 

1.8T0.7 

0162160501 

4.5T0.7 

1.8T0.4 

4.5T0.7 

1.8T0.4 


a Data are first modified by subtracting half of the minimum observed 
intensity for each direction and each line and then scaled for the observed 
variation of the HDFN flux and finally normalized by the model HDFN 
SWCX flux. 

b Scaled only by the observed variation of the HDFN flux. 

c Scaled by the observed variation of the HDFN flux and then scaled 
again by the model variation of the HDFN SWCX flux. 
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third observation to scale the data from the first 
and second observations, with the results listed 
in Table 2 and plotted in Figure 3a. However, 
our model prediction for the HDFN SWCX flux 
does vary by about 10% (see Table 5, and Fig- 
ure 6) with the intensity being lowest for the last 
observation. We use these model fluxes to add a 
second normalization to the SEP fluxes where the 
results are also listed in Table 2 and plotted in 
Figure 3a. There are a couple of caveats in us- 
ing this approach. First, we assume that the he- 
liospheric SWCX oxygen line emission dominates 
the observed flux. A significant contribution ei- 
ther from the Milky Way or yet more distant emis- 
sion or a constant flux from the heliopause would 
decrease the accuracy of the scaling by adding a 
constant offset. Second, we assume that the solar 
wind, both flux and abundances, on average are 
the same for both the northern and southern solar 
hemispheres. 

To gauge the affect of possible more distant 
emission we recalculated the normalized line in- 
tensities with the assumption that half of the min- 
imum observed fluxes in both directions and both 
spectral lines originates beyond the heliosphere. 
The results are listed in Table 3 and plotted in 
Figure 3b. Not surprisingly the sense of the varia- 
tion is the same as before and the only significant 
difference is that the fluxes are reduced by the fac- 
tor of two. 


Fig. 3. — a) The SEP line fluxes after using the 
HDFN data to normalize flux for the O VII (upper 
grouping of points) and O VIII line fluxes (lower 
grouping of points). The data are from Table 2. 
Data points with an X symbol are show the results 
when the HDFN SWCX flux is assumed to be con- 
stant. The points with a circle symbol include a 
further scaling to account for the model variation 
in the direction of the HDFN, see text for the de- 
tails. b) Same as a) except for the assumption that 
half of the minimum observed flux in both direc- 
tions (data from Table 3) has an origin outside of 
the heliosphere. 


2.3. Solar Wind Data 

The ACE and WIND satellites at the LI point 
provide density and velocity (among other quanti- 
ties) measurements of various species in the solar 
wind. Since SWCX emission is proportional to 
the flux of the solar wind species producing the 
. emission, the ACE and WIND data can be used 
to gain insight into the variation of the observed 
flux. Figure 4 shows the solar wind proton flux 
and the 0 +7 /0 +6 density ratio for the period of 
this program (the observation intervals are shown 
by the vertical lines). 

Both the solar wind proton flux and the oxy- 
gen density ratio show significant but not unusual 
variations over the interval. It is unfortunate that 
the data which are the most relevant, the actual 
0 +8 and 0 + ' fluxes, are not readily available from 
the A CE instruments and the two displayed para- 
meters must act as surrogates for at least a qual- 


7 





2003 Day of Year 


Fig. 4. — ACE solar wind proton flux (upper 
curve, units of 10 8 particles cm -2 s _1 ) and the 
0 +7 /0 +6 density ratio data (lower curve) covering 
the time interval of the observations. The vertical 
lines show the individual observation intervals. 

itative understanding of the circumstances. Note 
that both the solar wind flux and primarily the 
0+ 7 /0+ 6 density ratio are lower in the days lead- 
ing up to the third observation group than the 
similar periods before first and second observation 
groups. This supports the reasonableness for scal- 
ing the fitted O VII and O VIII fluxes. 

2.4. ROSAT All-Sky Survey Ecliptic Pole 
Data 

In addition to the XMM-Newton observations 
of the He focusing cone we also examined the 
RASS data from the ecliptic poles. Figure 5 shows 
the total keV count rate, the Long-Term En- 
hancement 2 (LTE, Snowden et al. 1997) and LTE- 
subtracted RASS \ keV surface brightnesses as a 
function of day of year for the south and north 
ecliptic poles (NEP). The data are from periods 
where the satellite pointing directions were within 
2° ofthe poles and are two-day averages. The re- 
moval of the LTE enhancements in essence applied 
a low-pass filter to the data with a time constant 
on the order of a day. While some of the tempo- 
ral variation in the '‘clean" data is possible due 

2 Long-Term Enalmcements are the name given to the obvi- 
ous SWCX enhancements which affected the RASS data. 
At the time (early 1990's) the origin of this temporally 
varying "background’’ enhancement was unknown. 



250 300 350 

1990 Day of Year 


Fig. 5. — Upper Panel : Two-day average RASS 
\ keV surface brightnesses for the intervals of the 
scans within 2° of the south ecliptic pole. The 
lower curve is the LTE surface brightness sub- 
tracted from the data (upper curve, Snowden et 
al. 1997). The middle curve is the “cleaned” 
(LTE subtracted) RASS surface brightness. Lower 
Panel: Same as the upper panel except the data 
are from the north ecliptic pole. Note that day 
366.0 of 1990 is day 1.0 of 1991. 





to variation in the cosmic background of distant 
origin (i.e., beyond the heliosphere), most is due 
to variations in the SWCX emission and the vari- 
able success in removing it from the RASS. Of 
note, however, are the enhancements in both the 
SEP and NEP data near day 345 in the total 
observed flux and a corresponding enhancement 
in the cleaned flux for the SEP. These are the 
strongest enhancements in terms of the integrated 
magnitude with the southern enhancement being 
somewhat stronger. Day 345 corresponds to the 
time period when Earth, and therefore ROSAT , 
were over the central direction of the He cone. The 
solar wind data in the early 1990s, the time period 
of the RASS observations (the last half of 1990), 
are too sparse (see Figure 9) to allow a signifi- 
cant quantitative analysis of the data. However, 
qualitatively the RASS data are in agreement with 
a model for enhanced SWCX emission during the 
period when ROSAT observed through the helium 
focusing cone. 

3. SWCX Model description 

The heliospheric SWCX model we use for our 
simulations is extensively described in Koutroumpa 
et al. (2006, 2007). This model is a self-consistent 
calculation of the solar wind charge-exchange X- 
ray line emission for any line of sight (LOS) 
through the heliosphere and for any observa- 
tion date, based on 3-dimensional grids of the 
IS neutral species (H and He) distributions in 
the heliosphere modulated by solar activity condi- 
tions (gravity, radiation pressure, and ionization 
processes). Highly charged heavy solar wind (SW) 
ions are propagated radially through these grids 
and the charge-transfer collision rates are calcu- 
lated for each of the ion species, including the 
evolution of their density due to charge-transfer 
with the IS atoms. With this process, we estab- 
lish 3-dimensional emissivity grids for each SW 
ion species, using photon emission yields com- 
puted by Kharchenko & Dalgarno (2000) for each 
spectral line following charge exchange with the 
corresponding neutral species (H and He individu- 
ally). Finally, the X-ray line emission is integrated 
along any LOS and observation geometry (for each 
observation date) in order to build the complete 
spectrum of SWCX emission in the given direc- 
tion. For comparison to present X-ray observa- 
tions we use the O VII triplet at 0.57 keV and the 


O VIII line at 0.65 keV, as they are the strongest 
spectral features and provide the best signal-to- 
noise ratio for the observations. 

We have conducted three simulations for each 
of the XMM observations that we name Models A, 
B and C. For all simulations the initial parameters 
of the IS neutral H and He flows are identical. The 
parameters specifying the IS neutral hydrogen are: 
n H (at 100 AU) - 0.1 cm" 3 , T = 13 000 K, V G = 
21 km s" 1 , A uw — 252.3°, (3uw — 8.5° (Lallement 
et al. 2005), where X uw and f3uw are the helio- 
ecliptic longitude and latitude respectively for the 
upwind direction of the incoming neutral H flow. 
The IS helium distribution has the following initial 
parameters: n// e (at 100 AU) = 0.015 cm -3 , T = 
6300 K, V 0 = 26.2 km s" 1 , \ uw = 254.7°, p uw 
= 5.3° (Witte 2004; Vallerga et al. 2004; Gloeckler 
et al. 2004). With those parameters, the He focus- 
ing cone is expected to be centered around (A dw, 
Pdw) = (74.7°, -5.3°). Model A assumes average 
solar minimum conditions (e.g., years 1996-1997 
or 2007-2008), Model B assumes solar conditions 
typical for late 2003, and Model C assumes av- 
erage solar maximum conditions (e.g., years 1990 
and 2000-2001). The basic differences in the three 
simulations are the neutral density distributions 
as they are shaped by the solar activity and ion- 
ization processes and the solar wind heavy ion dis- 
tribution. 

The H density reflects the action of gravity, ra- 
diation pressure, and losses due to CX collisions 
with SW protons and solar EUV ionization. The 
ratio, }jl , of radiation pressure to gravity for neutral 
H varies from 0.9 at solar minimum to 1.5 at solar 
maximum (Woods et al. 2000). The main source of 
ionization for H atoms is CX with solar wind pro- 
tons. Ionization rates are derived as a function of 
heliographic latitude from the SOHO/SWAN L y-a 
preliminary data analysis (Queinerais et al. 2006). 
For solar minimum the ionization rate decreases 
at around 40% from the solar equator to the solar 
poles, with an equatorial value of 6.7 x 10~ 7 s” 1 , 
while in solar maximum the relative decrease from 
the equator to the poles is only around 15% , with 
an equatorial value of 8.4 x 10" 7 s _1 . In late 
2003, with }i ~ 1.2. the equatorial ionization rate 
is 7.5 x 10~ 7 s" 1 and the equator to pole relative 
decrease is 30% in the H ionization rates. 

For He atoms the radiation pressure to gravity 
ratio is effectively zero and the atoms are grav- 
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itationally focused downwind (the helium focus- 
ing cone). The main cause of ionization is so- 
lar EUV radiation and electron impact. Recent 
work has shown strong evidence for a latitudi- 
nally anisotropic distribution of the He 30.4 nm so- 
lar irradiance and thus of the He photo-ionization 
rate (Witte 2004; Auchere et al. 2005) and it has 
also been shown that the electron impact ioniza- 
tion rate is also dependent on heliographic lati- 
tude (McMullin et al. 2004). In Koutroumpa et al. 
(2007) we tested preliminary grids of anisotropic 
He distributions in our calculation of heliospheric 
X-rays induced in SW T CX collisions. The present 
study, which examines emission from the He focus- 
ing cone, attempts to constrain the He ionization 
rates and their latitudinal dependance. 

For maximum solar activity, the photo-ionization 
rate is 1.4 x 10 -7 s -1 at the solar equator, while 
during solar minimum the photo-ionization rate 
at the solar equator is0.7xl0~ 7 s" 1 . The photo- 
ionization rates present a 50% decrease toward the 
solar poles. In 2003, the He photo-ionization rate 
was 1 x 10- 7 s -1 at the solar equator, again with 
a 50% decrease at the solar poles. 

The radial dependence of electron impact ion- 
ization is taken from Rucinski &; Fahr (1989) 
which is appropriate for solar minimum conditions 
but requires a three-fold increase for solar maxi- 
mum (Lallement et al. 2004). For this analysis we 
include a heliographic latitude dependence correc- 
tion factor, derived from McMullin et al. (2004), 
that we scale from minimum to maximum. The 
Lallement et al. (2004) scaling factor in combina- 
tion with the latitudinal anisotropy factor applied 
to the Rucinski & Fahr (1989) electron impact 
model is given in Table 4. An additional simu- 
lation (labeled Bl) was also performed which was 
identical to Model B except without electron im- 
pact ionization. 

The latitude dependence of the solar wind also 
affects the highly charged heavy ion distribution, 
where abundances depend on the solar wind type. 
During minimum solar activity, the solar wind is 
considered to be highly anisotropic, with a narrow 
equatorial zone (within ±20° of the solar equa- 
torial plane) of slow solar wind with an average 
speed of ~ 400 km s” 1 and the fast solar wind 
emitted from the polar coronal holes at a speed 
of ~ 700 km s -1 . The slow solar wind has a pro- 
ton density of ~ 6.5 cm -3 at 1 AU. while the fast 


flow is less dense at ~ 3.2 cm -3 at 1 AU. At so- 
lar maximum, the solar wind spatial distribution 
is considered to be a complex mix of slow and fast 
wind states that is in general approximated with 
an average slow wind flux. The ionic composi- 
tion of the two flows can be very different with 
the average oxygen content varying from [O /H] = 
1/1780 in the slow wind and [O/H] = 1/1550 in the 
fast flow. The charge-state distributions change 
as well, with the higher charge-states strongly de- 
pleted (or even completely absent, as for example 
O' 1-8 ) in the fast solar wind. For our model we 
adopt the oxygen relative abundances published 
in Schwadron & Cravens (2000): (0 +7 , 0 +8 ) = 
(0.2, 0.07) for the slow wind and (0 +7 , 0 +8 ) = 
(0.07, 0.0) for the fast wind, based on data from 
the Ulysses SWICS instrument (Schwadron et al., 
2000, unpublished document). 

The main difference in the SW heavy ion dis- 
tribution in our three simulations is the spatial 
(latitudinal) distribution of the slow and fast so- 
lar wind flows. For Model A (minimum SW) the 
slow SW is expanding in IP space through a ±20° 
equatorial zone on the solar surface, while the fast 
SW flow occupies the rest of the space. During so- 
lar maximum (Model C), we have considered that 
only an isotropic slow SW flow is present in the 
IP space. For the late 2003 period (Model B) we 
assume that there is no fast wind flow in IP space, 
in order to estimate average upper limits (we also 
assume that there is no CME or abnormal abun- 
dance distributions at the time of the XMM ob- 
servations) for the resulting SWCX X-ray emis- 
sion. Indeed, as demonstrated in Koutroumpa et 
al. (2006, 2007), for high ecliptic-latitude LOS, 
as is the case for the XMM He cone observations 
(/3 ~ -90°), the oxygen line intensity decreases 
from solar maximum to solar minimum conditions 
as the LOS crosses larger fast wind regions where 
the parent ions are strongly depleted. 

The model results for the SEP observations are 
listed in Table 5. In the same table we list also 
the HDFN simulation results for Model B. 

4. Discussion 

4.1. Comparison to XMM-Newton data 

In Figure 6 we compare the data and model re- 
sults for the oxygen line intensities for the three 
XMM-Newton observations. We can see that all 


10 


Table 4 

Heliographic Latitude Dependent Coefficients 


Hel. Lat. (°) 

Model A 
min (2008) 

Model B 
mid (2003-04) 

Model C 
max (2001) 

80 

0.6 

0.6 

2.4 

60 

0.65 

0.65 

2.4 

40 

0.8 

0.8 

2.7 

20 

0.9 

0.9 

2.85 

0 

1 . 

1 . 

3. 

-20 

0.9 

0.9 

2.85 

-40 

0.8 

0.8 

2.7 

-60 

0.65 

0.65 

2.4 

-80 

0.6 

0.6 

2.4 


N OTE . — Heliographic latitude dependent coefficients 
scaled for solar activity period to be applied on the Rucinski 
& Fahr (1989) electron impact ionization rate for He atoms. 


Table 5 

Model oxygen line intensities in LU 


ObsID 

Target 

Model A a 
0 VII 0 VIII 

Model B b 
O VII 0 VIII 

Model Bl c 
O VII O VIII 

Model C d 
O VII O VIII 

0162160101 

SEP 

1.7 

0.4 

3.3 

1.3 

3.4 

1.3 

2.5 

1.0 

0162160201 

HDFN 

- 

- 

3.2 

1.3 

_ 






0162160301 

SEP 

2.4 

0.7 

3.8 

1.5 

4.0 

1.5 

3.0 

1.2 

0162160401 

HDFN 

- 

— 

3.2 

1.3 







0162160501 

SEP 

1.9 

0.5 

3.6 

1.4 

3.8 

1.5 

2.8 

1.1 

0162160601 

HDFN 

- 

- 

2.9 

1.1 

- 

- 

- 




a Minimum solar wind conditions (e.g., 2008). 
b Medium solar wind conditions (e.g., 2003-04). 

c Medium solar wind conditions (e.g.. 2003-04) with no electron impact ionizations. 
d Maximum solar wind conditions (e.g., 2001). 
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Observation Observation 

Fig. 6.— a) Model and data line intensities for the O VII triplet. The data are from Table 2. Data points 
with an X symbol are show the results w'hen the HDFN SWCX flux is assumed to be constant. The points 
with a circle symbol include a further scaling to account for the model variation in the direction of the HDFN 
(double-triangle symbol). The model points are from Table 5. Inverse triangle symbols are for Model A, 
square symbols are for Model B, diamond symbols for Model B1 and triangle symbols are for Model C. 
b) Same as a) but for the O VIII line, c) and d) are the same as a) and b). respectively, except for the 
assumption that half of minimum observed flux in both directions has an origin outside of the heliosphere. 
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the simulations predict the maximum line emission 
on the second observation date, when the earth is 
located at 72° of ecliptic longitude. The nomi- 
nal position of the He cone is at ~75°, only three 
degrees from the second XMM-Newton observa- 
tion date. We also include the Model B prediction 
for the HDFN pointings during the observational 
campaign. If we take into account the HDFN vari- 
ation predicted in the SWCX model and scale the 
XMM-Newton data we find that the data yield 
the same trend as the model predictions, although 
with different relative line strength, both for the 
first scaling of the data and for where half of the 
minimum observed emission is assumed to be of 
cosmic origin. 

That the model absolute intensities underpre- 
dict the data (when no contribution from the cos- 
mic background is assumed) is not surprising since 
the observed field is not shadowed from distant 
emission, such as from the galactic halo. How- 
ever, assuming that the SWCX emission in the 
He cone is the only source of variation in the mea- 
sured data, we should obtain a residual Data (D) - 
Model ( M ) cosmic background (here after named 
R) that is constant over the time. This is not 
the case either for O VII or for O VIII. In order 
to have a constant background, that we assume 
to be equal to the average D - M difference be- 
tween the three observation dates for each simu- 
lation (R = | D — M |), we need to multiply each 
model point by a correction factor a: 

a = (D — R)/M, (1) 

which we show in Figure 7. Since, neither the aver- 
age R value representing the residual cosmic back- 
ground, nor the corrected SWCX emission (a ■ M) 
are allowed to have negative values, we assume 
that the average residual cosmic background can- 
not exceed the lower normalized value detected 
with XMM-Newton for each of the oxygen lines 
(7.2 and 1.7 LU respectively for O VII and O VIII ). 
A correction factor less than 1 (respectively more 
than 1) means that the model SWCX emission is 
overestimated (underestimated) and needs to be 
lowered (increased) in order to achieve a constant 
residual background. 

If we consider the second data-set. assuming 
that half of the observed flux is of cosmic origin, 
then the correction factor may be simplified as a' 
— D'/M. since D' is the residual observed SWCX 


emission, once the cosmic background has been 
removed. These correction factors are presented 
in the lower panels of Figure 7. 

From Figure 7 we find that the least correc- 
tion (best fit) is required for the simulations B/Bl 
(solar activity ajusted for late 2003) which is in 
good agreement with the data, especially for the 
O VII line. The model requiring the largest correc- 
tion is the one accounting for solar minimum ac- 
tivity period (model A), which assumes a highly 
anisotropic distribution for the SW oxygen ions. 
We also see that results for the O VII triplet need 
smaller corrections than the O VIII results, since 
the uncertainties both for XMM-Newton O VIII 
data and solar wind data for 0+ 8 are very large. 
Nevertheless, especially when assuming that the 
oxygen lines are dominated by the SWCX emis- 
sion from the heliosphere, the correction factor 
errorbars (calculated based on error propagation 
from the XMM-Newton data standard variations) 
overlap, revealing that all four models have little 
differences amongst them, with respect to the ab- 
solute differences from the data set. 

Our results suggest that toward high ecliptic 
latitudes in late 2003, there was little fast wind, 
and a slow wind flux with high oxygen content is 
required to reproduce the SWCX X-ray increase in 
the He cone. In order to investigate this assump- 
tion, we calculate the correction factors needed by 
model B to be scaled for different latitudinal ex- 
tents of the slow solar wind, and we plot the results 
for line O VII in Figure 8. 

Starting at a SW latitudinal extent of ±30° 
the correction factor progressively improves to its 
best fit for the three exposures, obtained for an 
isotropic slow wind (latitudinal extent of ±90° 
around the solar equator). However, any correc- 
tion factor variation remains within errorbars, and 
could eventually be considered non-significant. 

4.2. Comparison to ROSAT | keV data 

In Figure 9 we compare the ROSAT | keV sur- 
face brightness as a function of the observation 
date, with the equivalent SWCX model | keV sur- 
face brightness. The SWCX simulated spectra to- 
ward the North and South ecliptic poles, for each 
ROSAT observation date, were convolved with the 
ROSAT | keV band response in order to obtain 
the surface brightness due to SWCX emission for 


$ 


13 


Observation 


Observation 


Fig. 7. a) Correction factor a. to apply to the O VII model results in order to obtain a constant data-model 
difference corresponding to the residual cosmic background assumed this latter is constant, b) Same as a) 
except for O VIII model results, c) Correction factor a' to apply to the O VII model results in order to match 
the residual observed SWCX emission, assuming half of the minimum observed flux is of cosmic origin, d) 
Same as c) except for O VIII. The symbols have the same annotation as in Fig. 6. 





30 40 50 60 70 80 90 


Slow solar wind ± latitudinal extent (mean model) 

Fig. 8. — Correction factors to apply to model B results for different latitudinal slow SW extent. 

these dates and LOS. The first thing to note is 5. Conclusions 

that the model shows an expected enhancement 

for both the SEP and NEP directions. This is not We specifically designed a series of XMM- 

unreasonable as the He cone while centered along Newton observations in order to observe SWCX 

a line at -5.5° ecliptic latitude has an angular ex- emission from the helium focusing cone in the he- 

tent of ~ ± 25°. Comparing the model to the data liosphere. Three observations of the South Ecliptic 

shows a qualitative agreement with enhancements Pole, coupled with three monitoring observations 

in the total observed flux (upper curve), also in of the Hubble Deep Field-North were performed, 

both directions, and indeed they are the strongest The HDFN observations were used to normalize 

enhancements in terms of their magnitude and du- ^he SEP data, in order to attempt to eliminate 

ration. The cleaned data (lower curve) indicate short-term variation effects of the solar wind. We 

that the NEP enhancement was more effectively observe a SWCX heliospheric excess (with respect 

removed than that for the SEP 3 in the empirical °ff' cone directions) of 2.1 ±1.3 LU for O VII 

cleaning process. Considering both the cleaned anc ^ 2.0±0.9 LU for O VIII , when we assume that 

and uncleaned data, the SEP enhancement shows lines are dominated by SWCX emission from 

an asymmetry which is seen as a shoulder at higher ^he heliosphere. If we assume that half of the 

longitudes. This is a marginal result but is consis- observed emission is of cosmic origin, we find an 

tent with additional SWCX emission from the sec- SWCX excess of 1.7±1.2 LU and 1.2±0.7 LU for 

ondary He cone. However, this shoulder could also O VII and O VIII lines respectively. The maximum 

be the effect of a poorly timed enhancement of the emission is observed on the second XMM-Newton 

solar wind, as the He cone reacts fairly quickly to observation, when the observatory (earth) was lo- 

SW enhancements (see Figures 2 and 5 of Cravens cated at 72° of ecliptic longitude. This location 

et al. 2001), especially toward polar look direc- is consistent with the nominal central position of 

tions. the He cone (~74°). We find no clear evidence 

of the secondary He cone reported in Collier et 

3 Data from the southern hemisphere were more sparse than al. (2004), blit the data are not inconsistent with 

from the north as data losses from passage through the Such a component either. 

South Atlantic Anomaly preferentially affected the south, 

and this in turn affected the empirical cleaning process. COIlipaied the data results to a sclf- 

consistent model of the SWCX X-ray emission 
in the heliosphere, taking into account both H 
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Fig. 9.— Upper Panel: Model and data \ keV surface brightnesses for the south ecliptic pole. The plain 
curves are the total (upper) and LTE-corrected (lower) RASS ± keV fluxes. The diamonds are the SWCX 
model ^ keV fluxes for the same direction. On the right axis (dotted curve) we present dailv averages of the 
solar wind proton flux during the same period. Lower Panel: Same as the upper panel except the data and 
model results arc for the north ecliptic pole. Day 366.0 of 1990 is dav 1 of 1991. 





and He neutral components in the interplanetary 
space and also accounting for solar activity vari- 
ations. The model predicted an equivalent trend 
between the three observations, with the maxi- 
mum emission occuring at 72° of ecliptic longi- 
tude. The average absolute intensity predicted by 
the model is roughly half of the measured total 
intensity, since the field of view does not include 
any shadowing structures (i.e. dense molecular 
clouds) that would absorb the more distant halo 
emission. However, the on-cone excess predicted 
by the model (assuming solar conditions close to 
the late 2003 period, when the XMM-Newton ob- 
servations were performed) is much smaller with 
respect to the data on-cone excess, only 0.5 LU 
and 0.2 LU for O VII and O VIII respectively. 

This discrepancy between the model results and 
the data may be attributed to short-scale solar 
wind variations that the model was unable to ac- 
count for. The ACE and WIND data in the eclip- 
tic plane at 1 AU displayed significant variations 
during this period, but within reasonable limits. 
However, considering the high ecliptic latitude ob- 
servation geometry, both for the SEP and the 
HDFN pointings, we cannot assert that the solar 
wind conditions were indeed the same as the ones 
measured in the ecliptic, or that both the SEP and 
HDFN observations were influenced by solar wind 
of the same characteristics. In order to quantify 
the model to data discrepancy we produced a se- 
ries of scaling factors needed by the different sim- 
ulations to match the data. The corection factors 
for the simulations accounting for the late 2003 pe- 
riod of solar activity show a reasonable agreement 
of the model with the data. 

We also compared a series of RASS data to- 
ward the South and North ecliptic poles with the 
SWCX model results in the ~ keV band. The data 
(both NEP and SEP) exhibit a clear enhancement 
near the nominal He cone position predicted in 
the model. The SEP data also show a shoulder 
towards higher longitude which may be consistent 
with the presence of the secondary He cone. How- 
ever, it must be noted that the shoulder could also 
be due to a serendipitous enhancement of the solar 
wind. 

More X-ray observations of the He cone, espe- 
cially during different solar activity periods, are 
needed to further constrain the SWCX emission 
from the He cone and the heliosphere in general. 


Hopefully, the study will be reinforced by spectral 
diagnostics with the help of future X-ray missions 
including calorimeter instruments. 
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